use "C:\broz_data_ROIO_2011.dta", clear
*Model 4.1*
probit vote1 dwnom1  prezprty if  party1==0, rob cl(icpsr)
*Model 4.2*
probit vote1 dwnom1  prezprty if  party1==1, rob cl(icpsr)
*Model 4.3*
oprobit vote1 dwnom1  prezprty senate if  party1==0, rob cl(icpsr)
*Model 4.4*
probit vote1 dwnom1  prezprty senate if  party1==1, rob cl(icpsr)
*Model 4.5*
probit vote1 dwnom1  prezprty senate  h1945 h1959 h1962 h1965 h1976 h1978 h1980 h1983 h1998 s1959 s1978 s1980 s1983 s1991 s1998 s2009, rob cl(icpsr)
*Table 5 Substantive Results*
*First difference of DW-Nominate, Democrats*
estsimp probit vote1 dwnom1  prezprty senate if  party1==0, rob cl(icpsr)
setx mean
setx prezprty 0 senate 0
simqi fd(prval(1)) changex(dwnom1 -.2666664  -.0656328)
drop b1 b2 b3 b4
*First difference of DW-Nominate, Republicans*
estsimp probit vote1 dwnom1  prezprty senate if  party1==1, rob cl(icpsr)
setx mean
setx prezprty 0 senate 0
simqi fd(prval(1)) changex(dwnom1 .2769915 .4647024)
drop b1 b2 b3 b4
*First difference of Senate, Republicans*
estsimp probit vote1 dwnom1  prezprty senate if  party1==0, rob cl(icpsr)
setx mean
setx prezprty 0 senate 0
simqi fd(prval(1)) changex(senate 0 1)
drop b1 b2 b3 b4
*First difference of Senate, Democrats*
estsimp probit vote1 dwnom1  prezprty senate if  party1==1, rob cl(icpsr)
setx mean
setx prezprty 0 senate 0
simqi fd(prval(1)) changex(senate 0 1)
drop b1 b2 b3 b4
*First difference of President's Party, Republicans*
estsimp probit vote1 dwnom1  prezprty senate if  party1==0, rob cl(icpsr)
setx mean
setx prezprty 0 senate 0
simqi fd(prval(1)) changex(prezprty 0 1)
drop b1 b2 b3 b4
*First difference of President's Party, Democrats*
estsimp probit vote1 dwnom1  prezprty senate if  party1==1, rob cl(icpsr)
setx mean
setx prezprty 0 senate 0
simqi fd(prval(1)) changex(prezprty 0 1)
drop b1 b2 b3 b4
*Table 6 Natural Experiment*
*Model 6.1*
probit  vote1 dwnom1 prezprty senate, robust cl( icpsr)
*Model 6.2*
probit  vote1 dwnom1 prezprty senate atlarge senateXatlarge, robust cl( icpsr)
*Model 6.3 Senators Only*
probit  vote1 dwnom1 prezprty time_to_elec if senate==1, robust cl( icpsr)
*Table 7 Constituency and Special Interest*
*Model 7.1*
probit vote1 bankpac skills netm netx mktvalagprodii medinc h1983 h1998 if senate==0, robust cl(icpsr)
*Model 7.2*
probit vote1 dwnom1 prezprty bankpac skills netm netx mktvalagprodii medinc h1983 h1998 if senate==0, robust cl(icpsr)
*Model 7.3*
probit vote1  bankpac skills netm netx mktvalagprodii medinc s1998 s2009 if senate==1, robust cl(icpsr)
*Model 7.4*
probit vote1 dwnom1 prezprty bankpac skills netm netx mktvalagprodii medinc s1998 s2009 if senate==1, robust cl(icpsr)
*Model 7.5 Homogeneous States*
probit vote1 dwnom1 prezprty bankpac skills netm netx mktvalagprodii medin s1998 s2009 if senate==1 & homogenous==1, robust cl(icpsr)
*Table 8 Substantive Results*
*First difference of bankpac, House*
estsimp probit vote1 bankpac skills netm netx mktvalagprodii medinc h1983 h1998 if senate==0, robust cl(icpsr)
setx mean
simqi fd(prval(1)) changex(bankpac 2.375039 8.323187)
*First difference of skills, House*
setx mean
simqi fd(prval(1)) changex(skills .3141101 .4023941)
*First difference of net imports, House*
setx mean
simqi fd(prval(1)) changex(netm .1641294 0.2483097)
*First difference of agricultural production, House*
simqi fd(prval(1)) changex(mktvalagpr .0054122 0.0133625)
drop b1 b2 b3 b4 b5 b6 b7 b8 b9
*First difference of bankpac, Senate*
estsimp probit vote1 bankpac skills  netm netx mktvalagprodii medinc s1998 s2009 if senate==1, cl(icpsr)
setx mean
simqi fd(prval(1)) changex(bankpac 6.086522 15.794672)
*First difference of skills, Senate*
setx mean
simqi fd(prval(1)) changex(skills .2469205 .2861311)
setx mean
*First difference of net imports, Senate*
simqi fd(prval(1)) changex(netm .1263851 .1832335)
*First difference of net exports, Senate*
setx mean
simqi fd(prval(1)) changex(netx .0492512 0.0733254)
*First difference of agricultural production, Senate*
setx mean
simqi fd(prval(1)) changex(mktvalagpr .006889 0.014466)
drop b1 b2 b3 b4 b5 b6 b7 b8 b9
*end*
